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ABSTRACT 
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SIGNIFICANCE  AND  EXPLANATION 

Fluids  ar«  within  us  and  all  around  us,  from  kitchen  to  galaxy,  and  have 
been  studied  intensively  for  over  two  hundred  years.  Some  central  questions, 
however,  have  continued  to  defeat  Science.  Why  do  you  not  inhale  the  same 
body  of  stale  air  that  you  have  just  exhaled?  The  answer  lies  in  a  break-away 
of  fluid  streams  from  solid  surfaces  which  we  see  everywhere,  but  cannot  yet 
understand  or  predict  in  any  detail. 

It  now  appears,  however,  that  some  decisive  progress  towards  solving  the 
riddle  has  been  made  by  a  restricted  group  of  aerodynamicists  and  that  we  are 
in  the  midst  of  what  can  really  be  called  a  break-through.  .At  the  same  time, 
a  much  wider  scientific  community  appears  to  have  formed  a  strong  interest  in 
access  to  the  subject  on  account  of  its  importance  and  also,  of  dimly 
perceived,  novel  challenges  that  it  poses  for  mathematical  analysis. 
Conversely,  it  looks  as  if  mathematics  could  accelerate  and  widen  the  break¬ 
through.  To  this  end,  the  following  offers  an  introduction  to  the  main  line 
of  the  new  theory  for  a  wide  class  of  scientists  and  describes  some  of  its 
analytical  challenges. 

Since  the  earlier,  advanced  reviews  have  proven  so  difficult  for  all  but 
specialist  aerodynamicists,  the  following  is  restricted  to  just  the  basic 
theory  in  its  simplest  setting.  This  offers  an  opportunity,  also  here  seized, 
of  giving  a  clarified,  new  development  of  the  main  theory  from  very  few 
premises  to  illuminate  its  logical  structure. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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STEWART SON ' S  TRIPLE  DECK 


R.  E.  Meyer 

1 .  Introduction 

The  main  concern  of  boundary-layer  theory  during  its  first  half-century 
has  been  with  fluid  motion  at  large  Reynolds  number  Re  past  streamlined 
bodies.  In  this  context,  the  term  'streamlined*  is  meant  to  indicate  that  the 
motion,  except  very  close  to  the  body  and  in  a  thin  wake,  is  almost 
independent  of  Re  and  can  therefore  be  determined  closely  without  reference 
to  the  boundary- layer  development.  The  practical  scope  of  the  theory  is 
severely  limited,  however,  by  such  a  condition  of  'weak  interaction'  in  the 
sense  that  the  boundary  layer  and  wake  have  only  a  minor  influence  on  the  rest 
of  the  flow. 

The  primary  aim  was  to  understand  the  aerodynamics  of  wings  at  small 
incidence  (Figure  1)  and  observation  soon  confirmed  the  striking  success  of 
Prandtl's  theory  in  this  context.  As  the  incidence  is  increased,  however,  a 
stage  comes  where  the  airfoil  "stalls".  Indeed,  if  the  airfoil  of  Figure  1  be 
rotated  through  90°  to  a  position  broadside  to  the  incident  stream  (Figure  2), 
intuition  and  experience  predict  a  flow  pattern  in  which  the  stream  has 
"broken  away"  entirely  from  the  body  at  what  had  been  originally  its  nose  and 
tail.  At  an  intermediate  inclination  of  the  airfoil  (Figure  3),  the  incident 
stream  will  whet  the  front  part  of  the  upper  airfoil  side,  but  will  break  away 
from  it  at  some  point  which  inviscid  fluid  dynamics  by  itself  has  not  been 
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able  to  determine.  The  boundary  layer  on  the  body  and  the  rest  01'  the  flow 
then  'interact  strongly1  in  the  sense  that  no  approximation  to  either  is 
determinable  independently  of  the  other. 

It  is  readily  apparent  from  everyday  experience  that  a  similar  principle 
of  strong  interaction  applies  to  many  real  fluid  motions.  For  instance,  human 
breathing  must  involve  such  "breakaway"  of  a  stream  in  an  essential  way, 
because  otherwise,  we  could  only  inhale  the  same  body  of  gas  that  we  have  just 
exhaled.  Despite  all  the  aeronautical  triumphs,  therefore,  the  failure  of  the 
first  half-century  of  boundary- layer  theory  to  advance  beyond  the  weak- 
interaction  premise  has  left  the  understanding  and  prediction  of  most  real 
fluid  motions  still  in  a  rudimentary  stage. 

This  challenge  and  need  to  come  to  grips  with  strong  interactions  has 
fuelled  a  cooperative  research  effort  over  two  decades  under  the  leadership  of 
K.  E.  Stewartson,  which  has  resulted  in  a  major  step  forward.  The  following 
aims  to  open  a  readier  access  to  this  step  for  a  wider  audience,  because  one 
can  hardly  fail  to  be  impressed  by  the  work  achieved,  to  enjoy  the  elegance  of 
the  new  ideas  taking  shape,  and  to  be  intrigued  by  the  novel  challenges  to 
mathematics  emerging. 

The  present  survey  is,  of  course,  far  from  being  the  first  and  owes  a 
particular  debt  to  those  of  Messiter  [1978]  and  Stewartson  [1975].  The  early 
reviews,  however,  aimed  mainly  at  redirecting  research  and  were  addressed  to 
those  steeped  in  the  work,  in  the  first  place.  Later  reviews  were  burdened  by 
the  sheer  bulk  of  all  the  recent  results.  This  was  aggravated  by  the  fact 
that  the  earlier,  clear  successes  of  the  theory  concerned  supersonic 
aerodynamics,  where  the  theory  rests  on  a  rather  large  body  of  knowledge 
commanded  normally  only  by  professional  aerodynamicists.  To  avoid  this 
substantial  barrier  to  access,  the  introduction  to  follow  is  limited  to  low- 


spaed  flow,  with  only  a  few  remarks  on  the  supersonic  successes  in  Section 
8.  It  attempts  a  survey  taking  little  for  granted  beyond  the  Navier-Stokes 
equations  and  rudimentary  notions  of  fluid  dynamics  of  wide  familiarity.  And 
it  aims  only  at  an  illumination  of  the  main  features  of  the  rational  structure 
of  the  theory  to  the  degree  that  they  have  yet  emerged.  [Reference, 
accordingly,  is  made  only  to  books  and  review  articles,  rather  than  original 
sources.] 

The  next  section  outlines  the  limit  concept  underlying  Prandtl's 
boundary-layer  theory  for  weak  interaction.  Section  3  recalls  certain  aspects 
of  that  theory  which  are  of  particular  relevance  to  the  present  purposes. 

Those  begin  to  be  addressed  in  Section  4  by  the  introduction  of  a  triplet  of 
notions  which  appear  essentially  sufficient  to  characterise  the  new  theory. 
Section  5  shows  how  their  application  to  the  Navier-Stokes  equations  results 
in  a  need  for  more  than  two,  simultaneous  boundary- layer  limits.  A  third, 
needed  limit  is  developed  in  Section  6  and  shown  to  complete  a  definite 
triple-deck  of  limit  concepts,  but  all  the  same,  to  leave  us  far  short  of  a 
definite  system  of  differential  equations  for  a  description  of  the  flow. 
Section  7  outlines  how  an  exploration  of  the  interaction  between  the  limit 
concepts  ties  up  the  loose  ends  into  a  boundary  layer  problem  of 
mathematically  quite  novel  type.  It  has  been  explored  successfully  for 
supersonic  flow  (Section  8),  but  for  low-speed  flow,  turns  out  to  be  a  dead¬ 
end.  In  Section  9,  finally,  an  attempt  is  ma.de  to  sketch  Sychev's  ideas  on 
how  the  impasse  might  be  broken  to  nail  down  an  important  step  in  the 
understanding  of  fluid  motion. 

In  a  way,  the  whole  present  enterprise  may  be  described  as  an  attempt  to 
enable  the  non-specialist  reader  to  benefit  thoroughly  from  such  surveys  as 
Stewartson  [1975],  Messiter  [1978]  and  Stewartson  [1981].  To  that  end,  an 


Appendix  adds  a  succinct  derivation  of  the  Illingworth-Stewartson 
transformation  to  give  a  first  impression  of  how  that  transformation  removes 
the  barriers  to  a  direct  extension  of  triple-deck  theory  to  compressible  and 
supersonic  flow*  Its  developswnt  for  unsteady  and  three-dimensional  flow 
would  not  yet  appear  to  have  reached  a  stage  suitable  for  an  introductory 
review. 


2.  The  liinit  Concept 


To  fix  the  ideas,  the  discussion  will  be  focused  on  the  example  of  flow 
past  an  airfoil  such  like  that  indicated  in  Figure  1*  The  notion  will  also  be 
assumed  two-dimensional,  incompressible  and  steady  in  the  pilot's  frame  of 
reference,  to  which  it  will  be  referred  throughout.  Let  V  denote  the  flight 
speed  (the  wing  is  envisaged  as  moving  into  air  at  rest,  far  ahead)  and 
x,y,  Cartesian  coordinates  measured  in  units  of  a  body  dimension  L  from  a 
suitable  point  on  the  body  surface.  Mass  conservation  is  satisfied  by 
deriving  the  (non-dimensional)  velocity  components  u,v  in  the  respective 
directions  of  x  and  y  increasing  by 
(1)  u  ■  S'l'/dy,  v  -  -3^/3x 

from  a  stream  function  <|*(x,y)  measured  from  the  body  surface  in  units  of 

UL.  If  P,y  and  p(x,y)  denote  respectively,  the  constant  density  and 

2 

viscosity  and  the  local  pressure  measured  in  units  of  pu  ,  the  overall 
Reynolds  number  is  PUL/p  -  Re  and  momentum  conservation  is  expressed  by  the 
incompressible  Navier-Stokes  equations. 


♦ 
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y  xy 

x  yy 
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y 
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where  subscripts  denote  partial  differentiation,  with  boundary  condition 
u  -  v  «  0  on  the  body  surface. 

The  notion  of  limit  most  commonly  employed  in  Boundary-Layer  Theory,  and 
briefly  referred  to  as  Lim  hereafter,  is  that  of  a  'stretching'  transformation 

x  -  x/X(Re),  y  -  y/Y(Re) 

(3)  *(x,y)  -  o(Re)  <Mx*,y*) 

*  *  * 

p(x,y)  -  p#(x,y)  +  w (Re)  p  (x  ,y  ) 


-7- 


such  that  the  limit  Re  ♦  •  of  the  transformed  equations  with  appropriate 

•  *  *  *  *  * 

boundary  conditions  have  a  strict  solution  (x  ,y  ),  p  (x  ,y  ),  that  is, 

*  *  *  * 

p  is  differentiable  and  4*  is  thrice  differentiable  for  all  x  and  y 

in  the  fluid  domain. 

Superficially,  of  course,  this  notion  looks  quite  improper  in  the  sense 
of  routine  mathematics,  since  it  concerns  the  limit  of  the  differential 
equation  instead  of  the  limit  of  its  solution.  'Proper'  mathematical  work  on 
such  equations,  however,  has  long  been  proven  to  be  hopeless  until  an 
understanding  has  been  gained  of  the  real  issues  that  it  must  address. 
Meanwhile,  the  'improper'  limit  notion  has  proven  itself  a  thoroughly  reliable 
and  fruitful  guide  to  the  discovery  of  candidates  for  the  role  of  approximate 
solution  of  (2),  which  are  a  greatly  needed  starting  point  for  rigorous 
work.  Once  this  key  shift  of  objective  towards  the  discovery  of  candidates  be 
accepted,  the  theory  to  be  outlined  will  be  found  to  have  a  remarkable  degree 
of  logical  coherence,  even  though  the  present  account  is  not  arranged  to 
emphasize  that  aspect. 

One  obvious  example  of  such  a  Lim  is,  of  course,  the  'external  flow'  in 
which  x,  y,  ♦  and  p  themselves  are  kept  fixed  as  Re  ♦  “,  so  that  in  (3), 
(4)  X  =  Y  H  a  =  1,  ir  =  0 

and  ♦  *  ♦  (x,y),  p  ■  p  (x,y)  describe  the  classical,  irrotational  flow  past 
the  geometrical  airfoil  shape  with  circulation  determined  by  the  Kutta- 
Joukowsky  condition. 

Another  familiar  example  is  Prandtl's  boundary  layer,  for  which  it  is 
convenient  to  take  the  direction  of  x  increasing  tangential  to  the  airfoil 
surface  in  the  sense  of  the  local  flow  and  that  of  y  increasing,  normal  to 
that  surface  and  pointing  into  the  fluid.  Prandtl's  Lim  is  then  characterized 
by  the  stretching  transformation 
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of  the  local  coordinates. 

To  keep  confusion  between  different  Lira's  at  bay,  an  asterisk  will  refer 

to  the  generic  definition  (3),  but  a  subscript  i  will  denote  its  particular 

*  #  * 

case  (5),  e.g. ,  denotes  the  special  case  (5)  of  ♦  (x  ,y  ),  and 

v^  “  -3^/3x^,  etc.  The  subscript  e  only  serves  to  recall  the  shift  of 
emphasis  which  interpretation  of  (4)  as  a  Lin  introduces  into  the  aspect  of 
quantities  already  considered  in  (2). 

The  main  content  of  the  weak-interaction  principle  is  that  the  two  Lim's 
(4)  and  (5)  together  give  a  complete  picture  of  the  flow,  and  that  brings  up 
an  important  conceptual  point  intrinsic  to  the  limit  motion.  The  two  Lim's 
represent  different  views  of  the  same  solution  of  the  Navier-Stokes  equations 
(2),  different  magnifying  glasses  by  which  we  focus  on  different  features  of 
the  motion  so  as  to  gain  a  greatly  increased  understanding  of  what  goes  on  in 
the  fluid  by  the  comparison  and  correlation  of  the  two  views.  It  follows  that 
the  two  views  must  display  a  consistency  consonant  with  the  nature  of  the 
solution  of  (2)  which  they  illuminate.  The  diffusive  (parabolic)  structure  of 
(2),  in  turn,  implies  that  the  solution  of  (2)  must  be  a  smooth  function  of 
x  and  y. 

K  simple  illustration  of  such  consistency  is  provided  by  the  observation 
that  (1)  implies  v  (by  contrast  to  v^)  tends  to  zero. for  all  y^  under 
the  Prandtl  stretch  (5)  as  Re  -*■  ".  Since  v(x,y)  is  smooth,  its  external- 
flow  limit  v@(x,y)  must  also  be  zero  at  those  values  of  y  corresponding  to 
values  of  y^,  that  is,  at  y  =  0.  The  familiar,  zero-norraal-velocity 
boundary  condition  of  classical  hydrodynamics  then  re-emerges  as  a  consistency 
condition  implied  by  the  Lim-def inition  together  with  the  weak-interaction 


principle 


In  turn,  the  week-interaction  principle  is  also  meant  to  postulate  that 
the  external  flow  adheres  to  the  airfoil  without  breakaways  and  tends  to  a 
uniform  one  at  large  distances  from  the  body.  The  condition  of  zero  vft  at 
the  body  surface  then  suffices  to  determine  it  uniquely.  In  particular,  the 
values  of  p„  and  u.  along  the  body  surface  are  thereby  determined.  For 
the  pressure,  that  is  already  anticipated  by  the  notation  of  (3)  and  (4),  but 
for  the  tangential  velocity  component  u,  it  entails  a  further  consistency 
condition:  since  u^  *  u  under  the  Prandt 1-stretch  (5),  the  composite 
picture  of  the  two  Lim's  can  represent  a  smooth  solution  of  (2)  only  if 
u^  *  ue  for  all  x  at  those  y  belonging  to  both  Lim's,  i.e.,  as  *  m 
but  y#  “  y  0  (see  (6d)  below). 

Such  consistency  implications  of  the  limit  notion  play  a  key  part  in  the 
logical  structure  of  the  theory  to  be  outlined  below.  Their  technical  name  is 
matching  conditions,  and  their  mathematical  nature  has  been  explored  carefully 
by  Lagerstrom,  Kaplun  and  others  [Eckhaus  1979]  because  matters  can  get 
complicated  when  the  limit  notion  is  pushed  to  high  orders  of  approximation. 
For  the  present  purposes,  however,  no  such  complications  arise  and  it  may  be 
sufficient  to  mention  only  one  further  aspect:  if  consistency  of  two  Lim's 
turns  out  impossible,  then  there  are  only  two  alternatives  —  one  of  the  Lim's 
must  be  rejected  as  a  false  lead  on  the  way  to  discovery  of  a  candidate,  or  a 
third  Lim  must  intervene  between  the  first  two  to  reconcile  them. 
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3 •  Weak  Interaction 

A  fairly  full  account  may  be  found  in  many  books  [e.g. ,  Goldstein  1938, 
Lagerstrom  1964,  Stewartson  1964,  Meyer  1971]  of  how  the  limit  notion  a».d 
weak-interaction  principle  lead  uniquely  from  the  Navier— Stokes  equations  (2) 
to  a  set  of  limit-differential  equations  for  the  Lim  characterised  by  (5). 
Those  Prandtl  equations  are 

(6a)  u£  -  3'l»1/3yi#  vi  -  -34»i/3x1 

(6b)  ui3ui/3xi  +  vi3ui/3yi  -  -dp^/cb^  +  S^/Sy2 

(6c)  ui*xi'°*  -  vi(xi,0)  -  0 

(6d)  u1(xi,yi)  ♦  ua<xi>  **  yi  +  “ 

where  ue(x^)  is  the  surface  value  of  the  velocity  of  the  external  flow, 
related  to  its  surface  pressure  pe(x^)  by  Bernoulli's  equation 

(7)  pe  +  2  Ue  "  const* 

This  is  a  simpler  system  of  nonlinear  diffusion  equations  than  (2),  to  be 
complemented  by  specification  of  ui(0,yi)  for  all  yt  >  0  as  an  initial 
condition.  It  should  also  be  recalled  from  the  discussion  of  the  preceding 
Section  that,  from  the  viewpoint  of  the  Prandtl  Lim,  ue  and  pe  are  known 
quantities  in  (6b),  (6d). 

It  will  be  helpful  to  remark  here  on  three  further  aspects  of  prandtl 's 
Lim  which  are  of  particular  relevance  to  the  later  discussion.  While  our 
local  Cartesian  coordinates  have  a  strictly  global  interpretation  only  for  the 
limiting  case  of  an  airfoil  of  the  shape  of  a  flat  plate,  a  closer  exploration 
of  the  influence  of  surface  curvature  [e.g.  Goldstein  1938]  shows  it  not  to 
affect  the  theory  to  the  order  of  approximation  considered  in  this  account,  as 
long  as  the  body  shape  is  independent  of  Re  —  as  will  naturally  be  envisaged 
—  and  the  surface  curvature  is  finite  everywhere  except  at  a  trailing  edge. 
The  underlying  reason  flows,  of  course,  from  the  fact  that  a  Lim  like 


Prandtl's  uses  only  a  coordinate  system  covering  distances  y^  from  the  body 

surface  that  tend  to  zero  by  comparison  with  the  surface's  radius  of  curvature 

as  Re  ♦  **s  when  viewed  through  such  a  microscope,  the  body  surface  looks 

flat.  The  result  of  the  curvature  estimates  is  that  the  x^,  y^-coordinate 

system  can  serve  as  a  semi- local  one  for  all  the  Lira's  to  be  here  considered 

(except  the  external  flow)  in  which  x*  measures  distance  along  the  whole 

(upper  or  lower  side  of  the)  body  surface  and  y*,  distance  normal  to  it,  and 

that  interpretation  will  be  used  henceforth.  To  the  order  of  approximation 

considered  in  the  following,  moreover,  the  Lira-definition  can  be  applied  to 

(2)  without  regard  to  curvature  corrections. 

Secondly,  an  interesting  facet  of  the  weak  interaction  emerges  from  the 

consideration  of  consistency  to  a  second  approximation.  Though  v  -  -3<|i/dx 

tends  to  zero  under  the  Prandtl  Lim  (5)  as  Re  ♦  “,  the  same  does  not  follow 

1/2 

for  v^  «  -SiJij/Sx^  -  Re  v,  and  that  function  will  be  determined  whenever 

Prandtl's  equations  (6a)  -  (6d)  have  a  unique  solution.  Its  limit  v^ (x^,*) 

may  be  finite  and  non-zero,  in  which  case  "zero  normal  velocity"  ceases  to  be 

the  exact  boundary  condition  for  the  external  flow  at  finite  Re.  Instead, 

the  body  surface  assumes  a  'permeable'  aspect,  from  the  point  of  view  of  the 

external-flow  Lim  (4),  and  looks  as  if  there  were  a  small,  normal  velocity 
-1/2 

vfi(x,0)  -  Re  v^tx^,*)  across  it.  Its  determination  is  commonly  approached 
by  considerations  of  mass-flow  balance.  Since  (6c),  (6d)  show  u^  to  fall 
from  ue,  at  large  yi#  to  zero,  at  yi  -  0,  it  may  be  anticipated  that 
yiue*xi*  generally  exceed  the  stream  function  ^(x^y^,  although  both 

grow  with  yi  at  the  same  rate  as  yi  ♦  «,  by  (6d).  That  suggests 
consideration  of 

(8)  lim  (u  y  -  40/u  -  6(x.  ) 

v  -m»  e  ®  1 

yi 

y+0 
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which  characterizes  the  mass-flow  deficiency  of  the  boundary  layer,  and 
1  /2 

Re  6  ( x)  is  normally  found  to  be  a  well-defined,  positive  function  in  a 

weak  interaction.  In  its  terms,  the  surface-value  of  the  normal  velocity  for 
the  external  flow  becomes 

(9)  v  (x,0 )  -  (dfi/dx. )u  (x,0) 

e  i  e 

and  the  boundary-layer  mass-flow  deficiency  is  then  seen  to  be  interpretable, 

from  the  viewpoint  of  the  external-flow  Lim,  either  as  a  permeable  body 

surface  with  volume-outflow  rate  v_(x,0)  per  unit  distance  along  the  surface 

and  unit  span  normal  to  the  flow  plane,  or  as  a  solid-body  surface  displaced 

•1/2 

from  y  -  0  by  a  distance  Jlx^)  ■  0(Ra  ) .  Aerodynamicists  have  taken  to 

the  latter  interpretation  and  as  a  result,  mass-flow  effects  in  the  theory  are 

normally  discussed  in  terms  of  this  displacement  thickness  6. 

As  emphasized  in  the  preceding  Section,  the  known  arguments  leading  to 

Prandtl's  equations  (6)  prove  only  that  those  equations  identify  the  unique 

candidates  for  approximate  solutions  of  (2)  of  the  general  type  contemplated, 

if  there  are  any,  but  that  has  not  yet  been  established  rigorously. 

Historically,  the  value  of  Prandtl's  equations  was  first  demonstrated  by 

experiments  and  by  technology  based  on  them.  Applied  mathematicians  then 

explored  the  equations  quantitatively  and  purer  mathematicians  turned  to  the 

task  of  proving  that  their  solutions  are  approximate  solutions  of  the  Navier- 

Stokes  equations.  When  that  failed,  they  discovered  that  Prandtl's  equations 

themselves  have  mathematical  interest.  Oleinik  [1963]  proved  existence  and 

Halter  [1970]  found  a  substantially  simpler  proof.  Nichel  [1958]  established 

uniqueness  and  Serrin  [1967]  and  Pelletier  [1972]  attacked  the  key  question  of 

the  asymptotics  of  solutions  of  (6).  Pelletier  was  able  to  show  with 

considerable  generality  that,  in  normal  circumstances  and  for  dp  /dx.  <0, 

e  i 

the  influence  of  the  initial  condition  is  transient  and  with  increasing  xi# 
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the  solutions  develop  an  asymptotic  character  determined  by  the  functions 
p0(x^)  and  u0(x^)  related  by  (7).  In  this  way,  he  illuminated  the 
essential  character  of  Prandtl's  weak- interaction  equations  (6)  as  an 
inhomogeneous  system  of  partial  differential  equations  forced  by  the  given 
functions  dpe/dx^  and  u0(x^)  in  (6b)  and  (6d). 

All  those  theorems,  however,  amount  to  only  half  a  loaf  because  they 
depend  on  a  'favorable*  pressure  gradient,  dp^/dx^  <  0.  When  forced  by  an 
'adverse'  pressure  gradient,  dp0/dx^  >0,  the  solutions  develop  differently 
and  our  knowledge  is  more  tentative.  As  xi  increases,  the  surface  value  of 
9u^/dy^  then  decreases,  and  before  it  ceases  to  be  positive,  the  solution  of 
(6)  tends  to  develop  a  singularity  which  is  terminal  in  the  sense  that  a 
continuation  of  the  solution  beyond  it  is  not  possible.  Breakaway  (Figure  3) 
involves  a  reversal  of  flow  direction  near  the  airfoil  surface  and  hence,  a 
negative  'wall  shear'  du^/Sy^,  and  before  the  solution  of  (6)  can  reach 
that,  weak  interaction  appears  to  fail  decisively. 


4.  Mild  Interaction 


Strong  interaction  poaea  the  formidable  difficulty  of  interdependent  non- 
uniqueneae  of  both  the  external  inviscid  flow  (Figure  3)  and  the  viacoua  flow 
in  boundary  and  eh ear  layere.  It  ia  fortunate/  therefore/  that  careful 
experimenta  have  revealed  an  interval  of  'pre-atall'  incidencea  in  which 
recirculation  of  the  fluid  ia  confined  to  a  thin  'bubble'  region  adjacent  to 
the  airfoil  eurface  (Figure  4).  For  part  of  this  incidence  interval,  the 
obaerved  bubble  thickness  does  not  greatly  exceed  the  displacement  thickness 
of  the  weak-interaction  boundary  layer  ahead  of  the  bubble  and  accordingly, 
the  total  mass-flow  deficiency  near  the  upper  side  of  the  airfoil  remains  of 
the  same  order  as  in  a  weak  interaction,  and  the  first  approximation  to  the 
external  flow  therefore  remains  the  classical  potential  flow  past  the 
geometrical  airfoil. 

Such  observations  demonstrate  that  the  key-process  of  'separation'  — 
reversal  of  flow  direction  and  recirculation  of  fluid  close  to  the  body 
surface  —  can  occur  without  full  break-away  (Figure  3)  and  without  really 
strong  interaction  between  boundary  layer  and  external  flow.  An  opportunity 
thus  appears  of  dividing  the  difficulties  of  strong  interaction  by  studying 
first  the  flow  at  pre-stall  incidence  and  in  particular,  by  focusing  attention 
on  the  bubble  ends,  where  a  merely  local  breakdown  of  weak  interaction  seems 
to  occur.  This  opportunity  to  explore  the  local  mechanism  of  separation  was 
seized  by  K.  E.  Stewartson  and  the  aim  of  the  present  account  is  to  outline 
the  main  ideas  of  the  theory  so  far  developed  under  his  leadership. 

To  understand  a  little  better  how  he  succeeded  in  gaining  a  grip  on  the 
problem,  it  may  be  noted  that  the  bubble  thickness  (Figure  4)  is  not  observed 
to  vary  strongly  over  most  of  its  length.  This  indicates  that  the  rate  of 
change  of  displacement  thickness,  dfi/dx^,  also  remains  of  the  order  of 
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magnitude  typical  of  weak  interaction,  except  perhaps  at  a  bubble  end.  As 
expected  from  weak- interaction  experience,  on  the  other  hand,  confutation  of 
the  solutions  of  (6)  along  the  front  part  of  the  airfoil  surface  runs  into  the 
terminal  singularity  as  the  observed,  front  bubble  end  with  its  reversal  of 
flow  direction  is  approached.  That  would  seem  to  indicate  a  growth  of  > t- 
derivatives  of  velocity  or  pressure  to  levels  beyond  those  which  a  weak 
interaction  can  accommodate.  In  the  magnifying  glass  of  the  Prandtl  Lim,  this 
would  look  like  a  singularity,  and  its  resolution  might  require  a  stretching 
transformation  also  of  x.  These  plausible  interpretations  of  the 
experimental  and  theoretical  oracles  suggest  the  following  three  notions. 

(I)  The  displacement  thickness  remains  of  order  Re-1^2 
throughout. 

(II)  The  neighborhood  of  a  bubble  end  is  a  region  of  velocity  or 
pressure  changes  that  are  stronger  than  weak  interaction  can  describe, 
but  not  strong  enough  to  affect  the  external-flow  limit.  The 
possibility  needs  consideration  that  those  changes  might  penetrate 
farther  into  the  fluid  than  does  a  weak  interaction,  but  not  to 
distances  independent  of  Re  (which  the  external-flow  Lim  would  see). 

(Ill)  The  region  of  those  stronger  changes  near  a  bubble  end  appears 
quite  short  also  in  the  direction  along  the  surface,  by  comparison  to 
the  weak-interaction  boundary  layer  and  to  the  bubble  (Figure  4).  That 
exceptional  region,  in  which  the  apparent  singularity  is  resolved,  may 
plausibly  need  to  be  thought  of  as  shrinking  to  a  point  as  Re  ♦  “. 

One  of  the  aims  of  the  present  account  is  to  show  that  these  three 
notions  can  serve  as  the  essential  basis  of  the  theory  to  be  outlined  and  may 
therefore  be  considered  to  embody  a  definition  of  what  will,  for  brevity,  be 
called  mild  Interaction.  To  clarify  appeal  to  the  components  of  this 
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definition  in  the  following  discussion,  they  will  be  referred  to  briefly  as 
(I)  Mass-flow  bound, 

(ZZ)  Penetration, 

(ZZZ )  Localization. 

To  simplify  the  presentation,  it  will  be  very  helpful  to  restrict  its 
scope  by  three  working  principles  complementary  to  the  main  definition. 
Attention  will  be  confined  to  the  front  end  of  the  bubble,  which  is  easier  to 
think  about  because  it  involves  a  process  growing  out  of  a  standard,  un¬ 
separated,  weak- interaction  boundary  layer.  This  identifies  an  asymptotic, 
initial  condition  for  the  mild  interaction,  which  complements  the  definition 
and  will  be  referred  to  as 

(ZV)  Upstream  condition. 

Next,  the  notion  of  Lim  (Section  2)  admits  a  very  large  class  of 
stretching  transformations.  Zt  has  been  demonstrated  in  another  context 
[Meyer  and  Wilson  1973J  how  it  can  be  treated  systematically  to  sort  out 
uniquely  those  transformations  of  primary  significance  for  a  problem  at 
hand.  Such  an  approach,  however,  is  a  research  tool  for  settling  technical 
controversy  and  would  lengthen  the  present  discussion  so  as  to  obscure  the 
matters  of  real  interest.  As  a  working  principle,  an  a-priori  restriction  of 
the  stretching  functions  X,  Y,  0  and  if  of  Re  in  (3)  to  powers  of  Re 
will  therefore  be  adopted.  Zt  should  be  remarked  immediately  that  the 
advanced  theory  of  mild  interaction  bristles  with  logarithms  of  Re,  but  for 
the  present  purposes,  powers  will  turn  out  to  be  adequate. 

This  still  leaves  a  large  class  of  Lim's  and  not  all  of  them  can 
plausibly  be  of  equal  importance.  Many  will  have  to  be  discarded  as 
inconsistent  with  the  governing  equations  and  other  conditions  of  the  problem 
that  have  been  identified  already.  Others  will  be  merely  limits  of  other 


Lim's  that  offer  more  complex  and  useful  information  and  insight.  These 
aspects  of  the  Lim  notion  have  been  explored  by  Lagerstrom,  Kaplun  and  others 
[Eckhaus  1979]  and  it  would  again  be  unprofitable  for  present  purposes  to 
enter  upon  them  at  length.  Instead,  the  working  principle  will  be  adopted 
that  the  discussion  is  to  be  firmly  directed  towards  a  description  of  the  mild 
interaction  by  the  minimum  number  of  Lim's.  The  technical  term  for  them  is 
'significant  limits'  [Eckhaus  1979]. 
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5.  Two  Decks 

To  apply  now  the  notions  formulated  in  Sections  2  and  4,  begin  by  noting 
that  the  localization  condition  (III)  implies  a  stretching  function 

X  -  Re-<\  a  >  0 

in  (3)/  if  x  be  measured  from  a  suitable  point  in  the  mild  interaction 
region. 

The  upstream  condition  (IV)  requires  consideration  of  Prandtl's  Lim,  for 
-1/2 

which  Y  ■  0  ■  Re  ,  by  (5).  The  transformation  for  this  "main  deck"  Lim 
is  therefore 

x  *  xRe  ,  y  “  yli2 
(10)  ♦(x.y)  -  8^(x,y) 

p(x,y)  -  Pe(x)  +  i(Re)p(x,y) 

“1  /2 

with  9  -  Re  and  Pe(x),  u@( x)  used  as  in  Section  3  to  denote  the  known, 
surface  values  of  pressure  and  tangential  velocity  of  the  external  flow. 

Since  x  ♦  0  as  Re  ♦  •  at  all  x,  only  the  values  pe(0)  and  ue(0) 
actually  enter  into  the  argument  and  these  will  be  denoted  briefly  by  pe 
and  u@. 

Since  (10)  fails  to  account  for  the  penetration  condition  (II), 
consideration  must  also  be  given  to  another  transformation  that  can  bring 
greater  distances  from  the  body  surface  into  view,  and  this  is  called  the 
"upper  deck"  Lim.  At  such  greater  distances,  the  upstream  condition  (IV) 
implies  that  the  streamlines  came  from  the  external  flow  because  the  mass-flow 
bound  (I)  excludes  a  displacement  of  streamlines  near  the  body  surface  by  much 
more  than  $  -  0(Re  ).  By  Kelvin's  theorem,  therefore,  the  flow  in  the 
upper  deck  remains  an  irrotational  one,  which  has  no  distinguished  direction 
of  influence,  and  accordingly,  Lim  cannot  normally  involve  a  different  stretch 
for  x  and  y  and  must  have  a  transformation 


a  a 

x  “  xRe 


a  o 

y  “  yRe 


(11) 


♦<x,y)  -  u#y  +  a(Re)$(x,y) 


p(x,y)  ■  p  +  *<Re)p(x,y) 


with 

(12)  0  <  a  <  1/2 

in  order  that  the  upper  deck  do  penetrate  significantly  beyond  the  main  deck. 
The  condition  (I)  that  the  mild  interaction  has  a  mass-flow  deficiency 


u  y  -  +  of  the  same  order  as  that  of  a  weak  interaction  (section  3)  implies 
6 


(13) 


a 


-  Re~1/2 


Further  progress  depends  on  direct  use  of  the  Navier-Stokes  equations. 


If  D  abbreviates  the  convective  derivative  for  general,  stretched 
variables. 


(14) 


*3*3  3^*  3  3^i*  3 


+  v 


dx 


*  *  *  *  *  ' 

3y  3y  3x  3x  3y 


then  substitution  of  (10)  into  (2a)  yields,  by  (13), 


(15a) 


B  M  +  .  _Re-«p.  .  +  Re 

3y  3x  6  3y 


20-1  3  ^  j 

3x23y 


♦  0  as  Re  ♦  * 

by  the  Lim-definition  of  Section  2  and  by  (12).  Similarly,  substitution  of 
(10)  into  (2b)  and  of  (11)  into  (2a)  and  (2b)  yields 
(15b)  -B  3$/3x  +  if  Re^  20  3p/3y  +  0  , 

(16a) 

(16b) 

Mainly  through  the  localization  condition  (III),  the  enhanced  streamwise 
changes  in  the  mild  interaction  are  thus  seen  to  overshadow  the  viscous  shear 
in  both  decks  so  that  both  are  governed  by  inviscid  limit  equations  I 


A  A  A  A  1  —  2  O  .A  .A 

D  3i|i/3y  +  ir  Re  3p/3x  +  0 


-D  3<»/3x  +  t  Re1  20  3p/3y  ♦  0  as  Re  *• 
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To  learn  more  about  3  and  w  in  (15)  and  (16),  we  may  now  turn  to  a 
survey  of  alternatives  typical  of  Lim-arguments .  If  »  Re1  20  +  0  as 
Re  ♦  *,  (16a,b)  would  imply  mere  convection  of  a  velocity  perturbation 

already  present  upstream,  and  by  (IV),  could  describe  only  the  constant,  weak- 
interaction  correction  to  ue  and  vg  “  0 .  With  such  a  choice  of  stretching 
function,  therefore,  the  transformation  (11)  cannot  give  nontrivial 
information  on  the  velocity  perturbation  of  the  upper  deck.  If 
r  Re1  20  ♦  “,  on  the  other  hand,  (16a,b)  would  imply  p  =  const,  which  could 
describe  no  more  than  the  weak- interact ion  correction  of  pe  and  therefore, 
would  make  the  transformation  (11)  equally  uninformative.  The  only  stretch 
giving  nontrivial  information  can  therefore  be  n  =  Re2°  1 . 

Similar  information  on  tt  may  now  be  deduced  by  appeal  to  consistency. 

As  a  preliminary  remark,  recall  that  the  limit  equations  (16a,b)  are  an 
inviscid,  Euler  system  and,  by  the  upstream  condition  (IV)  and  Kelvin's 
theorem,  must  describe  a  potential-flow  perturbation.  This  could  be 
determined  uniquely  by  specification  of  the  'surface'  pressure  distribution, 
p(x,0),  together  with  the  condition  that  the  perturbation  decays  to  zero  as 
x2  +  $2  ♦  *.  Accordingly,  p(x,0)  f  0,  for  otherwise  the  perturbation  would 
vanish  identically  and  there  would  be  no  upper  deck  and  no  mild  interaction. 

By  (10)  and  (11),  however,  consistency  of  the  representation  of  pressure  as  a 
continuous  function  requires 

r  lim  $(x,y)  *  lim  p(x,y)  , 
y+0  y+* 

unless  another  Lim  intervenes  between  the  two  decks.  By  the  working  principle 
that  we  look  for  the  minimal  number  of  Lim's,  the  argument  against  such 
intervention  is  omitted,  since  p(x,0)  £  0  and  the  Lim-definition  requires 
p  to  be  bounded  on  compacts,  at  least,  it  follows  that 
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▼ 


(17)  *  -  n  -  Re2®-1  . 

By  (12)  the  limit  aquations  (15),  (16)  for  the  upper  and  main  decks, 
respectively,  now  become 

(18a)  D3$/3y  +  3p/3x  ♦  0  , 

(18b)  -&3ip/3x  +  3p/3y  +  0  , 

(19a)  03^/3y  ♦  0  , 

(19b)  -63^/3x  +  3p/3y  ♦  0 

as  Re  ♦  ",  for  all  x  and  x  and  for  all  y  and  y  >  0. 

The  plausible  suspicion  that  inviscid  equations  cannot  adequately 
describe  a  boundary-layer  process  can  be  raised  into  more  plastic  relief  by 
some  striking  conclusions  following  readily  from  (19a,b).  From  (19a), 

3$/3y  ■  u  is  merely  convected  through  the  main  deck  from  the  weak  interaction 
upstream,  so  that  it  is  a  function  only  of  ip, 

(20)  u  -  3ip/3y  -  f($)  . 

By  another  interpretation  of  (14),  it  is  also  seen  from  (19a)  that  the 
'streamline  slope',  defined  as 

(21)  -0*/3x)/(3*/3y)  -  9  , 

is  independent  of  y.  The  same  follows  for  30/dx,  and  by  (19b), 

30/3x  *  -(3ip/3y)  2  3p/3y  , 

so  that  3p/3y  is  also  Independent  of  y.  By  (10),  (11)  and  (17),  however, 

consistency  of  the  pressure  representation  requires 

lim  3p/3y  «  Re®  ^2  lim  3p/3y  +  0 

y+0 

by  (12).  Accordingly,  3p/3y  S  0  and 

(22)  0  and  p  are  independent  of  y  , 

the  main  deck  merely  transmits  the  perturbations  of  pressure  and  streamline 
slope  from  the  upper  deck  towards  the  body  surface. 


j 

i. 
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Now,  the  Navier-Stokes  equations  (2)  express  conservation  of  momentum  by 
a  balance  of  pressure,  momentum  and  viscous  stress.  By  their  boundary 
condition  at  the  body  surface,  the  momentum  near  it  must  be  arbitrarily  small, 
and  nontrivial  pressure  changes  in  a  mild  interaction  can,  by  (22),  be 
balanced  near  the  body  surface  only  by  another  deck  of  enhanced  shear  stress. 

Hindsight  now  permits  us  to  make  a  thumbnail  sketch  of  the  proof  of  this 
Section,  which  may  help  to  pull  its  threads  together  and  also  to  illuminate 
both  the  detailed  reasoning  behind  the  mi ld-interaction  premises  (I)  -  (III) 
and  the  possibilities  of  reducing  those  premises  further.  The  short  extent  of 
a  mild  interaction  prevents  viscous  shear  from  having  a  substantial,  direct 
influence  on  the  part  of  the  bubble-end  process  visible  in  the  classical 
Prandtl  limit.  As  a  result,  his  limit  describes  no  more  than  a  short  stretch 
of  weak  interaction,  unless  significant  perturbations  extend  beyond  his 
boundary- layer  thickness.  Nontrivial  perturbations  extending  beyond  it,  on 
the  other  hand,  make  the  process  visible  in  Prandtl' s  limit  one  of  mere 
transmission,  and  momentum  conservation  close  to  the  body  surface  is  not  then 
possible  without  enhanced,  viscous  shear  there. 


6 .  Lower  Deck 


Shear  stresses  substantially  exceeding  those  of  a  weak  interaction 
(Section  3)  require  Re  ^S^/Sy2  ♦  “  as  Re  ♦  •  and  the  third  deck 
exhibiting  them  must  be  a  Lim  with  transformation  of  form 


(23) 


x  -  x  Re°,  y  -  y  Re®,  3  >  ~ 

♦(x,y)  -  ff(Re)?(x,y),  p(x,y)  -  pg  +  if  p(x,y)  , 


with  again 

-  _  2 o-l 

x  -  Re 

for  consistency  with  the  transmitted  pressure.  Since  3  >  (10)  shows  that 

(23)  resolves  the  limit  y  ♦  0  of  the  main  deck. 

The  mass-flow  in  this  lower  deck  depends  on  the  particular  nature  of  the 
weak  interaction  from  which  the  mild  interaction  develops.  The  Reader  will 
have  grown  impatient  to  see  the  issues  resolved  in  one  case,  without 
enumeration  of  alternatives,  and  with  attention  already  focused  on  the  front 
end  of  the  bubble  (Figure  4),  it  will  now  be  further  restricted  to  the  most 
common  case  of  weak  interaction  with  normal  type  of  skin  friction, 
yi13<|»i/3yi  ♦  T^tx^)  as  T^(0)  ■  T,  briefly,  it  follows  from  the 

upstream  condition  (XV)  and  from  (10)  and  (12)  that  the  main  deck  must  also 
have 

(24)  lim  (y  *3^/3?)  ♦  t  as  x  ♦  -•  . 

y+0 


In  terms  of  $  and  u  =  3<P/3y,  this  relation  is 

lim  <if1/23*/3y>  *  (2t)1/2  , 

ip+0 


first  as  x  ♦  -®,  and  then  by  (20),  throughout  the  mild  interaction.  As 
y  ♦  0  and  y  ♦  ",  it  follows  from  (10)  and  (23)  by  consistency  that 
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a  Re^  3<>/3y  ~  (2t  a  i/8)1/2  as  y  ♦  " 
and  then  from  (13)  that 

(25)  3  *  Re  2®+1/2  an<j  3ij>/3y  -  Ty  +  A(x)  as  y  ♦  “  . 

Substitution  of  (23)  into  (2a,b)  now  yields,  by  a  computation  analogous 
to  that  of  Section  5, 


(26) 


-  3|  2 (a+0-1 )  3£  30-a-3/2  33j  ,  , 

D  “  +  Re  ~  -  Re  ♦  0 

3y  3x  3y 


3p/3y  *  0 


as  Re  +  Just  as  in  Section  5,  the  alternatives  may  now  be  explored  to  see 
fairly  readily  that  positive  or  negative  choices  of  a+3-1  or  3  -  1/2  -  a/3 
imply  truncations  of  (26)  that  are  either  degenerations  inconsistent  with  a 
mild  interaction  or  are  reduced  forms  of  (26)  describing  only  sublayers 
y  ♦  0  or  y  +  "  of  ( 23 ) .  The  choice 

o+B-1  -  3  -  j  -  o/3  *  0 

is  thus  seen  to  yield  the  only  significant  l.im  for  the  lower  deck,  and  it 
implies 

(27)  a  -  B  -  |,  a  -  Re~3/4,  i  -  f  -  «  -  Re" 1/4  , 

by  (25)  and  (17). 

All  three  stretching  transformations  (10),  (11)  and  (23)  have  thus  been 
determined  from  the  Navier-Stokes  equations,  the  definitions  of  the  limit 
notion  and  the  mild  interaction,  and  the  upstream  condition  (IV)  with  normal 
type  of  skin  friction.  Such  a  mild  interaction  must  penetrate  beyond  the 
normal  boundary  layer  thickness  of  order  Re-1/2  to  distances  of  order 
Re"3/8,  it  must  involve  three  significant  decks  and  must  exhibit  a  much 
enhanced  pressure  perturbation  of  order  Re-1/4  (by  contrast  to  that  of 
order  Re"1/2  in  a  weak  interaction). 
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From  (26)  and  (14),  the  lower  deck  must  be  governed  by  the  limit 
equations 

(28a)  u  -  3?/3y,  v  -  -3?/3x  , 

(28b)  u3u/3x  +  v3u/3y  ■  -3p/3x  +  3^u/3y^  , 

(28c)  3p/3y  -  0  , 

which  are  of  exactly  the  same  form  as  Prandtl's  equations  (6a),  (6b) l  Their 
theory  (Section  3)  leaves  no  doubt  that  the  physical  body-surface  condition  of 
the  Navier-Stokes  equations, 

(28d)  u(x,0)  -  v(x,0 )  ■  $(x,0)  -  0 

is  appropriate  to  the  lower  deck  equations  and  hence,  there  appears  no  call 
for  further  decks.  The  'outer  boundary  condition1  corresponding  to  (6d)  is 
now  (25),  and  since  the  lower  deck  in  a  weak  interaction  is  merely  the 
sublayer  close  to  the  body  surface,  the  initial  condition  for  (28)  is  (24) 
translated  by  (10)  and  (23)  into  its  lower-deck  representation, 

(28e)  u(x,y )  ♦  T  y  as  x  ♦  -«•  for  all  y  >  0  . 

Comparison  with  (25)  adds  the  information  that  A(x)  ♦  0  as  x  + 

With  demonstration  of  necessity  and  structure  of  the  triple  deck  thus 
complete,  it  may  be  opportune  to  vent  some  qualms.  Can  powers  as  small  as 
Re1/8  plausibly  distinguish  orders  of  magnitude  when  the  largest  Reynolds 
numbers  at  which  steady  motion  is  practically  realized  do  not  much  exceed 
IQ5?  My  own  view  is  that  the  true  asymptotic  structure  of  a  theory  tends  to 
emerge  only  at  a  late  stage  of  development.  It  is  reinforced  when 
experimental  results  at  even  lower  Re  show  agreement  [Stewartson  1981}  with 
theoretical  distinctions  based  on  even  smaller  powers  of  Re.  The  asymptotic 
notions  in  present  use  are  not  certain  to  be  definitive,  they  are  primarily  a 
fruitful  scheme  of  rational  guidance  to  the  discovery  of  candidates. 
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Another  qualm  Is  that  the  foregoing  proof  goes  to  show  that  the  candidate 
discovered  is  the  only  one  which  can  be  consistent  with  the  premises  outlined 
in  the  early  Sections.  It  does  not  address  the  question  whether  the  fairly 
complicated  system  of  equations  (18),  (19)  and  (28)  with  associated  conditions 
is  likely  to  be  consistent,  so  that  it  can  describe  an  actual  candidate.  In 
addition,  (28e)  sets  an  initial  condition  for  the  boundary- layer  type 
equations  (28)  which  is  anomalous  in  the  light  of  hydrodynamic  stability 
theory:  this  'boundary  layer'  is  to  start  out  with  a  'profile'  u(-*,y)  of 
which  every  point  is  an  inflection  point. 

A  more  urgent  concern,  however,  arises  from  uniqueness.  In  Prandtl's 
equation  (6b),  dp^dxj^  is  a  known  forcing  function,  but  no  information 
concerning  3p/3x  in  (28b)  is  yet  to  hand,  so  that  the  lower-deck  equations, 
at  this  point,  are  not  a  definite  set  of  differential  equations  at  all! 
Similarly,  ue(xA)  is  a  known  function  in  (6d),  but  the  outer  condition  (25) 
for  (28)  contains  an  unknown  function  A(x). 

The  upper-deck  equations  (18)  are  similarly  indefinite.  They  describe  a 
potential-flow  perturbation,  but  no  information  is  yet  to  hand  on  the  shape  of 
the  'body*  past  which  this  is  a  potential  flow.  Such  a  boundary  condition  for 
(18)  as  y  ♦  0  can  arise  only  from  the  main  deck,  but  that  has  been  seen 
(Section  5)  to  be  mainly  a  transmitting  mechanism,  so  that  the  missing 
condition  may  be  suspected  to  arise  actually  from  the  lower  deck.  The 
indeterminaciea  of  upper  and  lower  deck  may,  perhaps,  be  related? 


7 .  Interaction 


Apparently,  e  significant  interaction  between  the  decka  implicit  in  the 
mild-interaction  definition  or  limit  notion  remains  to  be  uncovered.  A 
natural  starting  point  is  the  transmission  statement  (22)  for  the  main  deck. 
By  (10),  (11)  and  (13),  the  upper-deck  representation  for  the  main-deck 


streamline  slope  is 


-8  - 


a?/a*  _ 

3$/3y 


*  a  -2°+7  A 
3^/3x  +  Re  u'y 

2  A  A 

u  +  Re  3$/3y 


-1  A  A 

♦  u^1 3$/3x 


as  Re  ♦  •,  by  (27).  Consistency  for  the  ratio  v/u  where  main  and  upper 
deck  merge  therefore  requires 

lim  3*/3x  ■  -u  lim  0(x,y) 

*  y~ 

By  (10)  and  (23),  similarly, 

-0  -  Re  *<3*/3x)/(3*/3y) 

and  some  information  on  this  ratio  is  obtainable  from  (25),  which  implies 

-1-2  -  - 

♦  -  2  T  y  "  A(x)y  ♦  B(x)  as  y  ♦  •  , 

whence 

(3*/3x)/(3*/3j)  A'(x)/x  as  y  ♦  -  , 

and  by  (27),  consistency  of  v/u  where  main  and  lower  deck  merge  requires 

0(x,O)  -  -Re"1/8A'(x)/T  . 

This  is  a  disappointing  result  because  the  transmission  statement  (22)  for 

streamline  slope  only  tells  us  now  that 

lim  3^»/3x  ♦  0  as  Re  ♦  • 
y+0 
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i.e.,  the  potent ial-f low  perturbation  of  the  upper  deck  is  a  flow  past  a  flat 
plate.  It  is  not  a  good  omen  for  the  search  for  a  non-vanishing  upper  deck. 


The  establishment  of  (27),  however,  has  changed  the  basis  for  the 
proofs.  While  the  limit  arguments  preceding  (27)  had  to  be  based  on 
comparison  of  powers  of  Re  restricted  only  by  inequalities,  an  opportunity 
now  arises  for  re-examining  the  limits  in  the  light  of  precise  knowledge  of 
those  powers  of  Re.  Indeed,  a  re-tracing  of  the  argument  leading  from  (15a) 
to  the  y- independence  of  (21)  —  verbatim,  except  for  the  addi tonal  fact  that 
a  -  3/8  --  shows  that  also 

5-  (r.,/8s>  -  -r.,/8(-21)"J  i3..  .s  r.  .  - 

3y  By  By 


1  /8 

so  that  Re  0  is  also  transmitted  unchanged  across  the  main  deck.  The 

preceding  calculation  now  supports  the  result 

(29)  lim  (Re1y/83$/3x)  -  u  A'(x)/T  , 

y+0  e 


which  expresses  a  mass-flow  interaction  between  the  lower  and  upper  decks. 

The  upper  deck  perturbation  is  therefore  a  potential  flow  past  a 

permeable  flat  plate  y  -  0  through  which  passes  an  outward  volume-flow  at  a 

rate  -Re  ^8u  A'(x)/T  per  unit  x-distance  (and  unit  span  normal  to  the  flow 
— —  e 

plane),  where  use  has  been  made  of  the  fact  that  x  *•  x  *  x,  by  (10),  (11) 
and  (23).  By  standard  potential  theory,  such  a  source  distribution  induces  on 
the  same  surface  y  »  0  an  x-component  of  velocity 


(30) 


»y 


u(x,0)  ■  -Re 


-1/8  ue  f-  A'(X1) 


—  / 

XT 


A 

x-x. 


dx 


1 


where  the  integral  denotes  the  Cauchy  principal  value,  if  it  exists  (and  x 
is  now  used  in  its  universal  meaning) . 


To  profit  frcm  the  other  half  of  the  transmission  statement  (22),  note 
that  (23)  and  (27)  cast  Bernoulli's  equation  for  the  potential  flow  in  the 
upper  deck  into  the  form 

P  +  j  (3^/3x)2  4-  j  (3i|»/3y)2  -  const 

-  P.  ♦  {  ~-,/8 *  1  (^)s  *  i  (^)21 

*  2  *  •  3$  2  3$  2  3$ 

4-  Re  1/2u'  y(3$/3x  4  Re  '^4u'  y) 


.  •  a* 

"  pe  +  2  ue 


which  is  its  value  as  x  ♦  This  is  again  of  no  direct  help,  but  a  re¬ 
examination  of  the  transformation  of  (2)  by  (11),  in  the  light  now  of  (27), 
shows  that  (16)  can  be  sharpened  to 


„  1/4  a 
Re  0 


3*/3y 

-3$/3x 


3p/3x 

3p/3y 


♦  0  as  Re  ♦  •  . 


If  the  classical  Bernoulli  calculation  is  applied  to  this,  it  supports  the 
stronger  statement 

Re1/,4(p  4-  j  (3i|»/3x)2  4-  j  (3<»/3y)2  -  pe  -  ^  u2]  -  0 
and  if  (23)  and  (27)  are  used  as  before  and  note  is  taken  of  (29)  and  (30), 


the  relation 


p(x,0)  ■  -Re1/®  u  u(x,0) 


of  form  familiar  from  linearized  potential-perturbation  theory  results.  This 
converts  (30)  into  an  equation  for  the  pressure  perturbation  on  the  'flat 
plate*  y  *  0,  and  since  the  pressure  stretch  is  the  same  for  all  decks,  by 
(27),  and  the  pressure  perturbation  is  transmitted  across  the  main  deck,  by 
(22),  it  is  also  an  interaction  condition  on  the  lower-deck  pressure. 
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So,  this  Hilbert  integral  equation  has  to  be  added  to  (28a-e)  and  (25)  to 


express  a  condition  of  mutual  consistency  between  the  three  decks. 

The  initial  similarity  between  Prandtl's  equations  (6)  and  the  mild-* 
interaction  equations  (28)  is  greatly  reduced  thereby,  and  the  relevance  of 
the  weak-interaction  knowledge  (Section  3)  for  triple-deck  theory  is  put  in 
doubt.  One  of  the  main  distinctions  is  that  Prandtl's  equations  (6)  are  an 
inhoangeneous  system  forced  by  a  known  external  flow,  but  the  mild-interaction 
equations  (25),  (28),  (31)  are  a  homogeneous  system  and  may  therefore  be 
suspected  of  admitting  non- trivial  solutions  only  for  particular  values  of  an 
as  yet  unidentified  eigen-parameter.  The  combination  of  boundary- layer 
equations  with  a  singular  integral  equation  poses  a  novel  mathematical 
challenge. 


8.  Successes  end  Failures 

It  will  be  recalled  that  the  objective  of  the  mild-interaction  concept 
was  to  resolve  local  breakdowns  of  the  weak- inter act ion  principle  and  most  of 
all,  to  elucidate  the  phenomenon  of  local  separation  characterized  by  reversal 
of  flow  directions  and  a  region  of  recirculation  of  the  fluid  close  to  the 
body  surface. 

One  of  the  simplest  breakdowns  of  weak  interaction  occurs  at  the  trailing 
edge  of  a  flat  plate  placed  edgewise  in  a  uniform  stream,  where  the  boundary 
condition  on  the  plane  y  •  0  changes  suddenly  from  no-slip  at  the  plate  to 
no-shear  at  the  wake-center.  Here  the  triple  deck  has  led  to  a  rather 
convincing  resolution  [Stewartson  1981]  of  the  weak-interaction  singularity. 
Admittedly,  it  is  not  nearly  as  simple  as  the  theory  here  discussed,  since  the 
analysis  appears  to  require  a  'staircase'  of  decks  [stewartson  1981],  but  the 
central  role  of  the  mild-interaction  concept  appears  well  assured.  For 
instance,  the  predicted  correction  to  Prandtl's  flat-plate  drag  fits  the 
experimental  observation  unbelievably  well  to  Reynolds  numbers  as  low  as  10. 
The  same  analysis  has  been  extended  to  wedge-shaped  trailing  edges  of 
sufficiently  small  wedge-angle  [Stewartson  1981].  While  those  are  welcome 
results,  however,  they  fall  short  of  addressing  the  main  issue  because  no 
separation  occurs  under  these  circumstances. 

A  much  more  resounding  success  has  been  achieved  in  the  analysis  of 
shock-boundary  layer  interaction  in  supersonic  flight.  Here  the  governing 
equations  are  much  more  complicated  than  (2)  because  the  fluid  motion  is 
intrinsically  coupled  to  thermodynamic  processes,  and  an  introduction  to  this 
branch  of  the  theory  would  completely  overload  the  present  account  with 
discussions  of  matters  not  central  to  the  illumination  of  the  triple-deck 
concept.  If  the  restrictions  be  accepted,  however,  which  support  the 
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Illingworth-Stewartson  transformation,  then  the  compressible  boundary-layer 
equations  for  weak  interaction  can  be  associated  with  a  system  of  Prandtl 
equations  (6)  in  a  way  explained  briefly  in  the  Appendix.  This  will  make  it 
easy  to  understand,  at  least  in  principle,  why  the  theory  set  out  in  the 
preceding  Sections  furnishes  an  almost  complete  guide  to  parallel  arguments 
extending  it  to  compressible-flow  problems.  What  emerges  [Stewartson  1975, 
Messiter  1978]  are  the  same  triple-deck  equations  (25),  (28)  and  this  will  be 
no  surprise  now,  for  they  govern  the  lower  deck  in  which  the  velocity  must  be 
small,  in  any  case.  There  is  one  drastic  change,  however:  the  supersonic, 
potential-flow  perturbation  in  the  upper  deck  is  governed  by  the  wave 
equation,  which  replaces  (30)  by  a  strictly  local  relation  leading  to  a  much 
simpler  interaction  condition, 

(32)  p(x)  =  -A  * ( x ) 

in  the  place  of  (31).  The  supersonic  triple-deck  equations  are  therefore 
still  a  homogeneous  system,  by  contrast  with  the  compressible  weak- interaction 
equations  (Appendix),  but  they  are  a  system  only  of  differential  equations. 

There  are  several  further,  favorable  circumstances  in  the  supersonic 
case,  and  the  simplest  example  for  illustrating  them  may  be  the  flow  past  a 
concave  corner  (Figure  5).  For  this,  the  wave  equation  and  indeed,  nonlinear 
inviscid  supersonic  theory  [Courant  and  Friedrichs  1948] ,  predict  a  very 
simple  solution:  two  regions  of  uniform  flow  separated  by  a  shock  springing 
from  the  very  corner  and  inclined  at  an  angle  related  directly  to  the  corner- 
angle.  Zn  reality,  this  can  only  be  the  asymptotic  pattern  far  from  the  body 
surface  because  the  effect  of  the  surface  boundary  condition  must  be 
transmitted  to  the  'external'  flow  by  a  boundary- layer  process,  on  which  the 
precise  location  of  the  shock  will  depend,  even  though  the  ultimate  shock- 
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inclination  does  not.  There  is  a  sense  in  which  this  is  a  strong  interaction, 
but  clearly,  the  adjustment  of  the  shock  location  is  relatively  easy  to  cope 
with. 

The  boundary- layer  process,  however,  poses  a  more  difficult  problem  when 
the  corner  is  sharp  (Figure  5).  For  if  the  radius  of  curvature  of  the  metal 
surface  is  much  smaller  than  the  displacement  thickness  (Section  3),  weak- 
interaction  theory  cannot  serve  and  mild-interaction  notions  are  the  least 
that  recourse  is  needed  to. 

Similar  considerations  apply  to  the  more  important  example  of  shock- 
boundary  layer  interaction  triggered  by  an  incident  shock  (Figure  6).  The 
external  flow  is  then  more  complicated,  but  can  still  be  synthesized  from 
simple  pieces  readily  adjustable  to  the  boundary- layer  process.  A  second, 
favorable  feature,  readily  apparent  also  in  the  corner-flow,  is  that  the  whole 
process  is  of  much  more  limited  extent  in  the  streamwise  direction  than  a  thin 
bubble  (Figure  4). 

The  most  important,  favorable  circumstance,  however,  was  pointed  out  at 
an  early  stage  by  Oswatitsch  and  Wieghardt  even  before  any  equations  relating 
to  such  interactions  had  ever  been  written  down.  The  paradox  was  that  the 
weak-interaction  equations  (Appendix)  are  parabolic  for  the  boundary  layer  and 
hyperbolic,  for  the  supersonic,  external  flow,  whence  mathematics  proves  that 
the  corner  (Figure  5)  can  have  no  influence  on  the  flow  upstream  of  it.  Just 
that  is  observed,  however,  and  to  help  explain  it,  they  pointed  out  that  the 
physical  influences  upon  each  other  of  the  mechanisms  of  supersonic  flow  and 
(weak-interaction)  boundary  layer  act  in  a  sense  favoring  a  self-contained, 
mutual  interaction  of  these  mechanisms.  In  first  translating  these  physical 
considerations  into  formulae,  Lighthill  was  able  to  give  strong  support  to  the 
plausibility  of  such  a  self-contained  process  indicating  an  eigensolution , 
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Stporotion  Dead  air  Attachment 

Fto.  I.  Sketch  of  main  properties  of  flow  when  a  shock  wave  interacts  with  a  laminar 
boundary  layer. 


6.  From  K. Stewartson,  Multistruct.ured  boundary  layers 
on  flat  plates  and  related  bodies,  Adv.Appl.Moch.  14 
p.  172,  Academic  Press,  New  York,  1974. 


several  properties  of  which  he  described.  This  gave  a  major  Impetus  to 


Stewartson's  formulation  of  the  triple-deck  concept  and  its  application  to 
shock-boundary  layer  interaction,  which  does  involve  separation. 

It  should  be  observed  that  the  attendant  reversal  of  flow  direction  near 
the  body  surface  introduces  a  striking  difficulty  so  far  ignored  in  the 
present  discussion.  The  boundary  layer  equations  (6)  or  (28)  are  parabolic 

with  a  nonlinear  diffusion  coeficient  related  primarily  to  the  tangential 

*  #  * 

velocity  component  u  =  /3y  (as  can  be  fully  clarified  by  the  Mises 

* 

transformation  given  in  the  Appendix).  In  a  weak  interaction,  u  >0 
everywhere  and  the  effective  diffusion  is  positive.  The  reversal  of  flow 
direction  inherent  in  separation,  however,  must  make  (28)  describe  a  nonlinear 
'heat-conduction*  process  that  turns  from  forward  to  backward  at  a  location 
which  is  obscure.  Clearly,  the  standard  coordinates  parallel  and  normal  to 
the  body  surface  are  then  unsuitable  variables,  but  discovery  of  a 
transformation  resolving  the  difficulty  at  an  affordable  price  does  not  seem 
imminent.  Meanwhile,  numerical  experiments  have  had  an  understandable 
tendency  to  degenerate  in*-o  a  long  brawl  with  unstable  diffusion.  .. 

In  addition,  shock- boundary  layer  interaction  is  a  much  more  complicated 
process  than  this  account  aims  to  discuss,  involving  at  least  two  triple  decks 
as  well  as  other,  distinct,  mild-interaction  processes.  All  the  same,  a 
combination  of  heuristic  asymptotics  with  numerical  experiments  has,  through 
the  efforts  of  a  number  of  investigators,  led  to  a  theoretical  prediction  of 
the  detailed  structure  of  supersonic  shock-boundary  layer  interaction 
(Stewartson  1975,  Messiter  1978] ,  which  displays  a  most  encouraging  measure  of 
agreement  with  experimental  observation.  Similar  success  has  been  achieved 
with  the  help  of  triple  decks  in  the  theoretical  description  of  supersonic 
flow  past  shallow  corners  (Figure  5)  and  even,  past  plate  sections  through 
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which  air  is  injected.  All  of  them  share  the  favorable  features  that  the 
process  is  self-contained,  that  the  external  flow  is  readily  mastered  and 
adjusted,  and  that  the  interaction  condition  is  (32).  To  learn  more  about  the 
successful  explanation  of  supersonic  flow  separation,  the  Reader  may  wish  to 
turn  to  [Messiter  1978],  tStewartson  1981]  and  the  literature  there  cited. 

The  subsonic  bubble-end  envisaged  at  the  start  of  the  present  account  as 
the  main  motivation  for  the  mild-interaction  definition  has  proven  much  more 
recalcitrant.  When  the  qualitative  influences  upon  each  other  of  the 
mechanisms  of  potential  flow  and  boundary  layer  are  considered,  a  change  of 
sign  is  found  to  occur  at  the  sonic  speed.  As  a  result,  the  physical  argument 
makes  a  self-contained,  mutual  interaction  implausible.  The  analytical  and 
numerical  exploration  failed  to  shake  this  prognosis  and  a  consensus  has 
developed  [Stewartson  1974]  that  the  triple-deck  equations  (25),  (28),  (31)  do 
not  possess  a  solution  describing  flow  separation. 


9.  Sychev's  Proposal 


To  break  this  impasse,  Sychev  suggested  that  some  of  the  mathematical 
foundations  of  the  theory  be  abandoned  by  relaxing  the  solution  concept 
(Section  2).  To  describe  his  proposal  better,  it  may  help  to  recall  first  an 
old  paradox  of  the  inviscid  break-away  theory.  To  fix  the  ideas,  flow  past  a 
cylinder  is  commonly  observed  (Figure  7)  to  break  away  so  as  to  leave  a  wide 
region  of  slowly  recirculating  fluid  as  a  wake.  The  first  of  all  the 
potential  flow  models  for  this,  due  to  Kirchhoff,  neglects  the  motion  in  the 
wake  so  that  the  pressure  there  is  a  constant,  pw>  That  reduces  the  non¬ 
uniqueness  of  the  potential  problem  to  that  of  a  family  of  definite  problems 
dependent  on  the  one  parameter  pw«  Wake  models  are  still  a  matter  of 
controversy,  but  Sychev's  considerations  involve  them  only  in  a  local  and 
conceptual  sense.  The  Kirchhoff  model  may  then  play  a  useful  role  as  a 
qualitative  representation  of  the  observation  that  the  pressure  on  the 
cylinder  surface  varies  much  less  rapidly  with  distance  downstream  from  the 
break-away  point  than  upstream. 

For  the  discussion  of  the  next  equation,  it  is  convenient  to  count  the 
arclength  x  along  the  body  surface  from  the  nose  stagnation  point  (Figure 
7).  The  separation  point,  denoted  by  x___,  marks  a  singularity  of  the 

Sep 

potential  flow  because  a  streamline  bifurcates  here;  the  branch  leaving  the 
body  surface  is  called  separation  streamline.  It  must  leave  the  surface 
tangentially,  because  a  stagnation  point  at  xsep  is  neither  plausible  nor 
observed,  and  potential  theory  shows  the  local  character  of  the  pressure 
distribution  on  the  surface  to  be  then  described  by  [Brodetsky  1923] 


(33) 


1/2 

p  ~  p  +  k(x  -  x)  as  x  t  x 

w  sep  sep 


w 


for  x  >  x 


sep 
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and  k  are 


with  constant  k.  Potential  theory  also  shows  that  pw,  xsep 
equivalent  parameters  standing  in  one-one  relation,  and  for  present 
purposes,  k  is  more  conveniently  regarded  as  the  family-parameter  of  the 
Kirchhoff  model.  Its  value  determines  the  potential  solution  uniquely  and 
k  >  0  is  found  to  imply  a  separation-streamline  curvature  at  xsep  exceeding 
the  body  surface  curvature  in  the  sense  that  the  separation  streamline  enters 
the  interior  of  the  body.  For  this  reason,  positive  values  of  k  must  be 
excluded  for  the  Kirchhoff  model. 

Negative  values,  on  the  other  hand,  would  imply  an  adverse  pressure 
gradient  of  arbitrarily  large  magnitude  upstream  of  *sepr  which  a  boundary 
layer  cannot  support  without  singularity  and  hence,  the  separation  must  occur 
before  its  defined  position  *8ep  is  reached.  This  contradiction  is 
sharpened  to  a  complete  dilemma  by  noting  that  k  =  0  would  imply  the  absence 
of  any  singularity  and  the  pressure  on  the  body  surface  then  continues  to  fall 
over  an  arclength  interval  extendng  beyond  xBept  so  that  the  weak- 
interaction  solution  also  continues  uniquely  and  there  is  no  break-away  of  the 
type  observed  (Figure  7).  This  paradox  inherent  in  the  Kirchhoff  model  was 
one  of  the  main  reasons  for  its  discreditation  over  many  years. 

Sychev  suggested  that  it  can  be  resolved  by  considering  k  to  depend  on 
the  Reynolds  number  in  such  a  way  that 
(34)  k(Re)  ♦  0  as  Re  ♦  •  . 

The  proposition  is  more  startling  them  may  be  immediately  evident.  It  shifts 
the  singularity  (33)  from  an  external-flow  property,  as  it  is  in  the  Kirchhoff 
model,  to  an  internal,  boundary- layer  feature.  An  attention-getting 
implication  is  that,  as  Re  +  00 ,  separation  is  predicted  to  take  place  at  a 
position  where  the  pressure  of  the  external  flow  is  still  falling  in  the 
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s^reamwise  direction.  That  le  not  only  quite  contrary  to  intuition,  but  also 
sharpens  the  negative  implications  for  subsonic  flow  of  the  physical  argument 
of  Oswatitsch  and  Wieghardt  (Section  8). 

The  relative  weakness  and  local  nature  of  the  singularity  (33),  (34) 
suggest  a  possible  connection  with  the  notion  of  mild  interaction.  The 
singularity  might  then  be  a  feature  of  the  potential-flow  upper  deck.  The 
next,  natural  question  for  Sychev  was  what  Reynolds-number  dependence  of  k 
would  be  consistent  with  this?  That  is  readily  deduced  from  (11)  and  (27), 
where  the  origin  of  local  coordinates  was  taken  at  an  appropriate  point  in  the 
mild-interaction  region,  now  seen  to  be  xsep.  Accordingly,  x  -  xsep  is  to 
be  interpreted  as  just  x  and  pe,  identified  with  pw.  In  the  upper-deck 
notation  (11),  the  relation  (33)  then  becomes 


and  by  (27),  since  the  respective  stretches  for  x  and  p-p  are  the  same  in 

© 


all  three  decks. 


k  - Y  Re”1/16  as  Re  ♦  * 


p - Y  |x| 1/2  as  x  ♦  -■ 


as  x  ♦ 


with  constant  Y  >  0,  because  positive  values  of  k  are  excluded  for  the 
same  reason  as  in  the  Kirchhoff  paradox. 

It  may  be  worth  remarking  that  (33)  and  (35)  together  express  a  condition 
of  consistency  between  the  triple-deck  Lims  and  the  ever-present,  fourth  deck 
of  the  external-flow  Lira.  What  the  observer  looking  at  the  external  flow  sees 
as  a  local  limit  (33)  describing  a  singularity  (at  finite  Re),  must  be  seen 


by  the  observer  looking  at  the  upper  deck  as  an  asymptotic  statement,  but  not 
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neccessarily,  as  a  more  detailed  condition.  In  Sychev's  proposal, 
accordingly,  the  role  of  the  triple  deck  is  to  resolve  the  apparent 
singularity  of  the  external-flow  Lim  needed  to  make  the  total  picture 
consistent  with  a  locally  Kirchhoff  model  of  break-away. 

However,  (35)  contradicts  (31)  because  it  makes  that  integral  roundly 
divergent.  Sychev  proposes  to  avoid  that  difficulty  simply  by  an 
interpretation  of  the  integral  in  the  finite-part  sense.  The  more  explicit 
meaning  of  (31)  then  becomes 


(36) 


p(x)  +  Y  u2|x|1/2H(-x) 


%  A'  (s)+YT  sV2H(s) 
7TT  •*— •  - 


X-S 

where  H(x)  denotes  the  Heaviside  unit  step  and  the  integral  is  again  in  the 
sense  of  Cauchy's  principal  value. 

The  introduction  of  the  finite  part  involves,  of  course,  a  substantial 
shift  of  the  solution  concept  into  the  realm  of  distributions  and  thereby 
mandates  a  revision  from  the  ground  up  of  the  theory  outlined  in  this 
account.  A  rational  motivation  of  the  Sychev  form  of  mild-interaction  theory 
however,  is  not  nearly  as  urgent  a  concern  as  the  look  forward,  and  in  a 
striking  role  reversal,  it  was  now  Stewartson  [1974)  who  called  publicly  for 
mathematical  existence  proof  1  To  appreciate  the  substantial  and  novel 
character  of  this  challenge,  it  may  help  to  collect  here  the  set  of  equations 


(25),  (28), 

(35),  (36)  arrived  at: 

(37a) 

u  *  3i/3y, 

v  ■  -3^/3x  , 

(37b) 

u3u/3x  +  v3u/3y  = 

-dp/dx  +  32u/3y2  , 

for  all  x 

and  all  y  >  0  , 

(37c) 

u(x,0)  »  v(x,0 )  - 

$(x,0)  for  all  x 

(37d) 

u  -  Ty  ♦  A(x)  as 

y  ♦  *,  for  all  x 

with  known  constant  T 


(37e) 


for  all  y  >  0 


u  ♦  Ty  as  x  ♦ 

(37*)  |x|_1/2  p  ♦  -Y  as  x  ♦  -• 

with  unknown  constant  Y  >  0, 

(37g)  p  ♦  0  as  x  ♦  " 

and 

(37h)  p(x)  +  Y  u2|x| 1/2H(-x)  - 

•  x-s 

In  the  sense  of  Cauchy's  principal  value,  with  H(x)  denoting  the  Heaviside 
unit  step.  It  is  still  a  homogeneous  system  for  the  unknown  functions 
u(x,y),  v(x,y),  p(x),  A(x),  and  the  uppermost  question  is  now  clearly  whether 
an  eigenvalue  Y  exists  for  which  these  equations  possess  a  nontrivial 
solution? 

Sychev  was  content  with  showing  that  such  a  solution,  if  there  is  one, 
would  have  asymptotic  properties,  as  x  ♦  •,  that  give  a  plausible 
description  of  the  flow  reversal  and  recirculation  inherent  in  separation. 

The  existence  question  was  addressed  by  Smith  [1977],  and  in  the  absence  of  a 
mathematical  theory  for  anything  resembling  (38a-h) ,  he  attacked  it 
numerically.  It  is  not  hard  to  appreciate  that  the  task  was  very  difficult 
and  could  not  be  accomplished  with  definitive  completeness,  but  Smith  [1977] 
did  succeed  in  establishing  very  substantial  evidence  in  favor  of  the 
conjecture  that  an  eigenvalue  of  Y  exists  near  0.44. 


Appendix:  Xllingworth-Stewartson  Transformation 
Interdisciplinary  research  is  made  appreciably  more  difficult  by  the 
custom  that  the  literature  of  each  branch  of  science  takes  a  great  deal  for 
granted.  Thus  reviews  of  triple-deck  theory  pass  quickly  to  the  supersonic 
case  where  the  main,  early  successes  are  found  (Section  8).  In  so  doing,  they 
take  a  remarkable  correlation  between  compressible  and  incompressible  boundary 
layers  for  granted,  and  it  may  help  non-aerodynamicists  to  gain  ready  access 
to  the  more  advanced  reviews  of  triple-deck  theory,  if  a  succinct  exposition 
of  that  correlation  be  here  appended. 


(i)  Boundary  Layer  Equations 

The  Havier-Stokes  equations  for  steady,  two-dimensional  flow  are 
3(Pu)/3x  +  3(Pv)/3y  <*  0 


x  3x  3y 

3v  .  3v.  3p  3  3u  .  3vx1  .23  r  3v  3ux1 

p(u  37  +  v  ^  "  ’  if  +  37  lV%  +  37)]  +  3  3?  [W(2  37  ' 

_  ,  3t  3t.  3p  3p  .  3  3t.  3  ,  >  3t. 

-  r_,3u,2  .  _,3v,2  ,3u  .  3v.2  2  ,3u  ,  3v.2, 

*-y[2(^)  +  2  ("5y )  +<37+37)  "3(3^+37)] 

together  with  the  equation  of  state  of  a  perfect  gas 

-  p  *  (c  -c  )PT  . 

P  v 

Here  cp,  cy  denote  the  specific  heats,  assumed  constant,  X,  the  heat 
conduction  coefficient,  and  T,  the  temperature.  For  a  fairly  full 
background  of  these  equations,  reference  may  be  made,  e.g. ,  to  [Meyer  1971) 
Only  weak  interaction  will  be  considered  in  this  Appendix.  For  a 
boundary  layer  on  a  solid  surface  y  ■  0,  a  direct  extension  to  these 
equations  of  the  transformation  (3),  (5)  reduces  them  to 


(38a) 


3(Pu)/3x  +  3(pv)/3y  »  0 


(38b) 

(38c) 


...  3u  .  3u.  dp  3  ,  3u. 

p(u  51  ♦  v  1J1  ’  '  &  +  a?  Iu  1J> 

„  .  it,  3t.  do,  I  ,  it.  ,  ,3u,2 

PVU  e"!;1  "r?1*!;""!;1 


(38d)  (c  -c  )PT  -  p(x)  ■  p  (x) 

p  v  e 

as  Re  *  U  L  Pg/Ug  ♦  ",  where  the  asterisks  and  subscripts  of  Sections  2,  3 

have  now  been  omitted  and  U,  L,  PQ  and  pQ  denote  reference  values  of 

velocity,  length,  density  and  viscosity,  respectively,  and  the  variables  are 

nondimensional .  The  standard  boundary  conditions  at  a  solid  surface  are 

(38e)  u  *  v  -  0  and  T  or  3T/3y  specified  at  y  *  0  for  all  x 

and  the  conditions  of  consistency  with  the  external,  inviscid  flow  are 

(38f )  u  ♦  u  (x),  T  ♦  T  (x),  p(x)  ■  P_(x)  as  y  ♦  "  , 

0  6  6 

where  ue,  T#,  pe  denote  again  the  body-surface  values  of  the  external 
flow.  These  are  Prandtl's  compressible  boundary- layer  equations. 

As  in  the  low-speed  case  (Section  3),  a  finite.  Re-independent  curvature 
of  the  body  surface  does  not  affect  the  equations  to  the  order  in  Re 
considered,  and  x  and  y  in  (38a-f)  can  therefore  again  be  thought  of  as 
curvilinear,  orthogonal  coordinates  along,  and  normal  to,  the  body  surface. 


(ii)  Mises  Transformation 

Let 

Pu  *  3i(i/3y,  pv  *  34>/3x,  ♦(x,0)  *  0 

to  satisfy  (38a)  and  take  as  new,  independent  variables  the  stream  function 
<Kx,y)  and  the  distance  s(x,y)  along  streamlines.  In  the  Lim  under  the 
stretch  (3),  (5), 

s(x,y)  -  x 

and  for  any  differentiable  function  f(^,s),  the  first  derivatives  become 


3f 


3f  „  3f  .  3f  3f 
3^  "  -pv  3*  +  17  '  37“  PU  3? 


Accordingly,  the  momentum  equation  (38b)  is  transformed  into 
(39) 

This  form  of  the  local  momentum  balance  raises  the  roles  of  convection  and 
nonlinear  diffusion  into  clearer  relief  and  is  often  useful  in  technical  work 


3u  dp  3  ,  3u. 

PU  77  “  "  dJ  +  PU  H  iVPU  3?}  ' 


(iii)  Crocco  Relation 

In  the  absence  of  chemical  reactions, the  stagnation  enthalpy, 

1  2  2 

c  T  +  —  (u+v),  of  inviscid  flow  is  constant  along  streamlines  and  in 
P  2 

particular,  therefore,  its  body-surface  value  in  the  external  flow, 

1  2 

c  T  +  —  u  *  const.  *  I  . 
p  e  2  e  e 

For  gas-dynamical  background,  reference  may  be  made,  e.g. ,  to  [Meyer  1971). 
1  2 

If  I  -  CpT  +  —  u  denotes  the  local  value  of  the  same  quantity,  then  the 
usual  combination  of  the  energy  equation  (38c)  plus  u  times  the  momentum 
equation  (38b)  reads 


(40) 


3  , .  3t  .  3u. 

3?  (X  37  +  uu  37> 


3  .  3l 

37  (P  37} 


if  the  Prandtl  number  uc^/X  -0*1.  This  assumption,  made  in  the  following 
is  a  common  gas  dynamical  approximation  because  the  range  of  o  is  normally 
very  small  and  not  too  far  from  unity;  for  air,  e.g.,  0  =  0.72, 

approximately,  over  quite  a  wide  range  of  thermodynamic  conditions. 

If  there  is  also  no  heat  transfer  to  the  body  surface,  the  boundary 
condition  there  is  3T/3y  -  0,  and  by  (38e),  also  3l/3y  =  0  there.  An 
obvious  solution  of  (40)  satisfying  this  boundary  condition  and  consistent 
with  (38f )  is 

I  =  const.  -  I 

e 
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Assuming  uniqueness ,  this  establishes  a  simple,  direct  dependence  of  local 
temperature  on  local  velocity, 

(41)  cT  -  I  -  4  u2  , 

p  e  2 

which  eliminate  the  need  for  further  consideration  of  (38c). 


(iv)  Xllingworth-Stewartson  Transformation 

If  U  *  T,  which  also  adopted  now  as  a  common  gas  dynamical 

0.89 

approximation  (for  air,  e.g. ,  u  «  T  over  quite  a  wide  range  of 

conditions),  then  by  (38d),  UP  is  also  constant  across  the  boundary  layer, 

i.e.,  UP  depends  on  x,  but  not  on  This  clearly  simplifies  (30)  a 

little,  and  if  account  is  also  taken  of  the  Bernoulli  relation,  p  u  du  /ds  « 

e  e  e 

-dp/d s,  for  the  inviscid  external  flow,  (39)  takes  the  form 

3u  PeUe  due  3_  .  jhx. 

3s  pu  ds  PePe  3<»  u  3t|» 


If  now  u  is 

then  3u/3a  - 

(42a) 


referrd  to  the  local  unit  u.  by  setting 

© 

u(<>,s)  -  u  (s)  u‘ (<t>,s)  , 


u  du'/ds  +  u'du  /ds  and  the  momentum  equation  becomes 
6  6 


3u' 

3s 


,.2- 


u'u 


p  /p  du 
e  e 

ds 


3  ,  ,  3u' . 

We  3?  (U  dT) 


with  boundary  conditions 

(42b)  u' (♦.s)  ♦  1  as  ♦  ♦  •  for  all  s  , 

(42c)  u'(0,s)  ■  0  at  the  body  surface,  for  all  s 

This  is  a  system  involving  only  u'(<P,s)  and  functions  of  s  because 

p  /p  •  t/*T  ,  by  ( 38d ) ,  and  on  account  of  the  Crocco  relation  (41). 

e  e 

Mow  recall  the  equations  (6)  for  the  incompressible  boundary  layer,  which 
are  seen  to  be  the  special  case  of  (38a,b,e,f)  resulting  for  constant  p  and 
M#  aay»  P^  and  u^  respectively.  If  now  the  Mises  transformation  to 
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independent  variables  is  applied  to  (6b)  and  if  *  ueiui'  then 

the  resulting  equation  is  the  special  case  of  (42)  for  p  »  p@  =*  p^,  p^  =  p^: 


(43a) 


3u.  u’2-1  du  .  a  3u' 

i  ,  i  ei  3  ,  ,  _ i, 

3s  u!u  .  ds .  ^iPiUei  3t|>  ui  3i(». 
i  l  ei  i  i  l 


(43b) 


u[  -  1 


as  i|>  ♦  “  for  all  , 


(43c) 


u|(0,si)  =0  at  the  body  surface,  for  all  s^ 


Now  (42a-c)  and  (43a-c)  are  seen  to  be  the  same  equation  system,  if  we 


identify 


♦  *  ♦i  , 


(44)  p  P  u  ds  -  p.p.u  ,ds.  =  d£,  say, 

e  e  e  i  i  ei  1 

2 

u'  -P  /P  .  du  ,  du  . 

. .  e  -1  e  -1  ei 

(45)  - - -  u  — r  -  u  .  — -  . 

,2  ,  e  d(  ei  ds. 
u*  -1  i 

That  (45),  like  (44),  refers  only  to  the  external  flows  can  be  deduced  from 
Crocco's  relation  (41)  as  follows.  By  (38d), 


„  1  2 
V  2  U 


„  A  1  2  1  ,2  2 

Ve  2  ue“  2  u  Ue 


Ue  2 

1  +2Tt-  (1-U’  > 

p  e 


2  2  2  2 

so  that  u*  -  p  / P  *  u'  -  1  -  (1-u'  )u  /(2c  T  )  and  (45)  becomes 
e  e  p  e 


u  .  du  «  du  . 

r .  ^  e  ■>  -1  e  -1  ei 

2c  T  Ue  Ai  *  Uei  d£  * 
p  e 


This  completes  proof  of  the  Theorem  discovered  almost  simultaneously  by 
Illingworth  and  Stewartson: 

Theorem.  To  any  inviscid,  compressible  external  flow  without  chemical 
reaction  and  with  surface  velocity  distribution  ue(x),  we  can  associate  an 


inviscid,  incompressible,  external  flow  of  surface  velocity  distribution 
uei(xi)  by 
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(47) 


duei/dxi 

U  P  u  ^ 
i  ei 


u  du  /dx 

+  — - — )  — S - 

2c  T  ’  2 

p  e  UP  u 
e  e  e 


and  for  0  =  1,  u  *  T  and  no  heat  transfer,  the  weak-interaction  boundary 
layers  are  then  related  by 

(48)  u!  -  - f-  -  u’  -  Hiliii  . 

1  Uei  Xi  Ue  x 

The  theorem  furnishes  the  solution  of  the  compressible  problem  as  soon  as 
the  solution  of  the  associated  incompressible  problem  has  been  obtained, 
though  first  in  the  Mises  form  u  =  u(ip,x).  To  complete  the  quasi-Cartesian 
form  of  the  solution,  a  quadrature  suffices: 

Corollary  1.  At  fixed,  corresponding  x  and  x^, 

dy  i  i 

dy.  P  UT 
i  e  e 

Proof.  At  fixed  s  and  s by  definition,  d<>  =  pudy  *  Piu^dy^  so  that 
dy/dy^  ■  P1u±/(Pu)  and  the  Corollary  follows  from  (38d). 

It  is  worth  remarking  that  the  theorem  can  be  simplified  by  reference  to 


the  local  Mach  number 


Me(x)  *  ue(x)/ae(x) 


where  a  (x)  is  the  speed  of  sound  defined  by 


(49)  a  »  Yp  /P  =  (Y-1)c  T 

e  e  e  p  e 

with  now  Y  =  c  /c  ,  by  (38d).  By  (41)  and  (49), 
p  v 

cpdTe/dx  =  _uedue/dx  ' 


dM  /dx  du  /dx  dT  /dx 
e  e  e 


2c  T 
p  e 


du  /dx 
)  — 5 - 


so  that  (46)  becomes  simply 


Uei/Me  =  const<  “  a0 


say,  and  (44)  completes  the  proof  of 
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Corollary  2.  (47)  in  the  theorem  may  be  replaced  by 


uei(xi>  "  a0  Me(x)  ' 


‘i  '  (WiPia0r1  f*  Ue{s)Pe 


( s )  a  ( s )  ds 


Corollary  1  may  also  be  simplified  to 


Vi 

a  P 

e  e 


<»i  *  ¥  ¥ 'o1  <’  -  ui2,ari1 


with  integral  taken  at  constant  x  because,  at  corresponding  points  (♦/X^) 
and  (<i*,x),  the  theorem  gives  u^/u  -  uei/ue  ”  a0^ae'  by  (50),  and  T  is 
given  in  terms  of  u  by  (41). 
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